A CONSTRUCTIVE BOREL-CANTELLI LEMMA. 
CONSTRUCTING ORBITS WITH REQUIRED STATISTICAL 

PROPERTIES. 

STEFANO GALATOLO, MATHIEU HOYRUP, AND CRISTOBAL ROJAS 

Abstract. In the general context of computable metric spaces and com- 
putable measures we prove a kind of constructive Borel-Cantelli lemma: given 
a sequence (constructive in some way) of sets Ai with effectively summable 
measures, there are computable points which are not contained in infinitely 
many Ai. 

As a consequence of this we obtain the existence of computable points which 
follow the typical statistical behavior of a dynamical system (they satisfy the 
Birkhoff theorem) for a large class of systems, having computable invariant 
measure and a certain "logarithmic" speed of convergence of Birkhoff averages 
over Lipshitz observables. This is applied to uniformly hyperbolic systems, 
piecewise expanding maps, systems on the interval with an indifferent fixed 
point and it directly implies the existence of computable numbers which are 
normal with respect to any base. 
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1. Introduction 

Many results in mathematics ensure the existence of points satisfying a given 
property P by estimating the measure of P and proving that it is positive. In 
general this approach is not constructive and does not give an effective way to 
construct points satisfying the given property. 

A key lemma in this kind of techniques is the well-known Borel-Cantelli one: 
Borel-Cantelli Lemma. Let {An} be a sequence of subsets in a probability space 
{X,fi). If fi{A„) < oo, then /i(limsupA„) = 0, that is, the set of points which 
are contained in infinitely many A„ has zero measure. 

Under these conditions, X — limsup A„ is a full measure set and hence it contains 
"many" points of X. In this paper we give a general method to construct points 
in this set. This method will be applied to some nontrivial problems, as construct- 
ing numbers which are normal in every base and typical trajectories of dynamical 
systems. 

To face this problem we will put ourself in the framework of computable metric 
spaces. Let us introduce and motivate this concept. It is well known that the state 
of a physical system can be known only up to some finite precision (because of 
measuring errors, thermal shaking, quantum phenomena, long range interactions 
etc.). From a mathematical point of view this knowledge is represented by a ball 
with positive radius in the metric space of all possible configurations of the system. 

In practice, the knowledge of the state of the system up to some finite precision 
can be described by a sentence like "the position of the point in the phase space at 
time 3 is a;(3) = 0.322 ± 0.001". What is important here is that it admits a finite 
description (a finite string of characters). 

This finite string of characters, can then be elaborated to estimate, for example 
the position or the distance of the system's status at time 3 with respect to other 
points of the space. 

This kind of identification 

Strings <-> [Points, Geometrical objects] 

if often implicit, and considered to be obvious but it underlies the concept of Com- 
putable Metric Space. 

A Computable Metric Space is a metric space where a dense countable set (which 
will be called the set of ideal points) is identified with a set of finite strings, in a 
way that the distance between points in this set can be computed up to any given 
approximation by an algorithm having the corresponding strings as an input (see 
section 12.21 for precise definitions) . 

For example in M the set Q can be identified with the strings "p over q" in a way 
that the distance between rationals can be obviously calculated by an algorithm 
having the strings as input. We remark that if K is considered as a computable 
metric space, then beyond Q there are many other points which admit finite de- 
scriptions, for example tt or are not rationals but they can be approximated 
at any given precision by an algorithm, hence in some sense this points too can 
be identified to finite strings: tt for example can be identified with the finite pro- 
gram which approximates it by rationals at any given precision. This set of points 
is called the set of computable real numbers (they were introduced by Turing in 
|Tur36| ) . The concept of computable point can be easily generalized to any com- 
putable metric space. Coming back to our main question, now the problem we face 
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is the following: Given some property P about points of X (or equivalently a subset 
of X), can this property be observed with a computer? That is, docs there exist 
computable points satisfying this property? 

For instance, given a (non atomic) probability measure let P be a subset of 
X of probability one: a point chosen "at random" will almost surely belong to P. 
But, as the set of computable points have null measure (is a countable set) the full 
measure of P induces a priori nothing upon its computable part (i.e. the set of 
computable points belonging to P). 

We will give some results which give a positive answer to this question when P 
is constructed by a Borel-Cantelli technique. Let us illustrate this (for a precise 
statement see theorem [1]): 

Theorem A. Let us consider a sequence of closed sets {An)nefi (with some effec- 
tivity condition, see definition^ such that ^/i(A„) < oo in an effective way (see 
Def U^. 

If the measure fi is computable (Def. ^) then there are computable points outside 
limsup that is lying in An 's only finitely many times. 

Computable absolutely normal numbers. As an example, a classical question 
where this kind of tool can be naturally applied is the normality: given a fixed 
enumeration base b of real numbers it it quite easy to prove that the set of b- 
normal numbers (the numbers where all the digits {0, 6—1} appear with the same 
frequency) has Lebesgue-measure one. Can we find computable normal numbers? 
The construction proposed by Champernowne [Cha33j happens to be algorithmic, 
so it gives a positive answer to the question. 

A natural and much more difficult problem is to construct numbers which are 
normal in every base (see sec. l4.2l for some historical comments on the problem). In 
section the existence of computable absolutely normal numbers will be obtained 
as a quite simple corollary of Theorem A. 

Computable points having typical statistical behavior. The above result 
on normal numbers is a particular case of the construction of computable points 
which follows the typical statistical behavior of a dynamical system. We will need 
the notion of computable dynamical systems, let us introduce it. 

The notion of algorithm and computable function can be extended to functions 
between computable metric spaces (Def. IB]). This allows to consider computable 
dynamical systems over metric spaces (systems whose dynamics is generated by 
the iteration of a computable function), and computable observables. With these 
definitions, all systems which can be effectively described (and used in simulations) 
are computable. 

Computable points (as described above) are a very small invariant set, compared 
to the whole space. By this reason, a computable point rarely can be expected to 
behave as a typical point of the space and give rise to a typical statistical behavior 
of the dynamics. Here, "typical" behavior means a behavior which is attained for a 
full measure set of initial conditions. Nevertheless computable points are the only 
points we can use when we perform a simulation or some explicit computation on 
a computer. 
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A number of theoretical questions arise naturally from all these facts. Due to 
the importance of the general forecasting-simulation problem these questions also 
have a practical importance. 

Problem 1. Since simulations can only start with computable initial conditions, 
given some typical statistical behavior of a dynamical system, is there some com- 
putable initial condition realizing this behavior? how to choose such points? 

Such points could be called pseudorandom points. Meaningful simulations, show- 
ing typical behaviors of the dynamics can be performed if computable, pseudoran- 
dom initial conditions exislQ. A somewhat similar problem has already been inves- 
tigated in |KST94j in the setting of symbolic dynamics. They consider recursive 
discretisations of the system (that is a subset of computable points) and look for 
conditions to ensure that a finite observer is unable to distinguish the motion on 
the recursive discretisation from the original system. 

In our framework, a first topological result is the following: if the system is 
computable and has at least a dense orbit, then there is a computable point having 
a dense orbit (see Thm. |3|). 

From the statistical point of view we can use the above Theorem A to prove the 
following second main result which we summarize informally below (see Thm.s [2] 
and [5] for precise statements). 

Theorem B. // {X,ii,T) is a computable dynamical system and 

(1) PL is a computable invariant ergodic measure. 

(2) The system {X,T,p) is log^ ergodic (see definition ] 15\ ) for observables in 
some suitable functions space. 

Then there exist computable points x for which it holds: 



for any continuous function / : X — > M with compact support. 

The above theorem states that in such systems there are computable points 
whose time average equals the space average for any such observable on X , hence 
providing a set of computable points which from the statistical point of view behave 
as the typical points of [X, fi) in the Birkhoff pointwise ergodic theorem. 

We remark that the approach taken in ^KST94j is quite different, in the sense 
that they give sufficient conditions (in terms of Kolmogorov complexity) for a subset 
of computable points (a recursive discretisation) which ensure that this set satisfies 
a kind of finite ergodic theorem (a much weaker property than II. ip but give no 

^It is widely believed that computer simulations produce correct ergodic behaviour. The 
evidence is mostly heuristic. Most arguments are based on the various "shadowing" results(see 
e.g. |HK95I chapter 18). In this kind of approach (different from our), it is possible to prove that 
in a suitable system, any "pseudo" -trajectory, as the ones which are obtained in simulations with 
some computation error is near to a real trajectory of the system. 

So we know that what we see in a simulation is near to some real trajectory (even if we do 
not know if the trajectory is typical in some sense). The main limit of this approach is however 
that shadowing results hold only in particular systems, having some uniform hyperbolicity, while 
many physically interesting systems are not like this. 

We recall that in our approach we consider real trajectories instead of "pseudo" ones and we 
ask if there is some computable point which behaves as a typical point of the space. 
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method to construct such computable points (because these conditions cannot be 
verified in a constructive way). 

To apply theorem B to concrete systems the main difficulty is to verify the 
points 1) and 2). In section [5] we show that these are verified for the SRB invariant 
measure (the natural invariant measure to be considered in this cases) in some 
classes of interesting systems as uniformly hyperbolic systems, piecewise expanding 
maps and interval maps with an indifferent fixed point. 

2. COMPUTABILITY 

The starting point of recursion theory was to give a mathematical definition 
making precise the intuitive notions of algorithmic or effective procedure on sym- 
bolic objects. Every mathematician has a more or less clear intuition of what can 
be computed by algorithms: the multiplication of natural numbers, the formal 
derivation of polynomials are simple examples. 

Several very different formalizations have been independently proposed (by Church, 
Kleene, Turing, Post, Markov...) in the 30's, and have proved to be equivalent: they 
compute the same functions from N to N. This class of functions is now called the 
class of recursive functions. As an algorithm is allowed to run forever on an in- 
put, these functions may be partial, i.e not defined everywhere. The domain of a 
recursive function is the set of inputs on which the algorithm eventually halts. A 
recursive function whose domain is N is said to be total. 

We now recall an important concept from recursion theory. A set i? C N is called 
recursively enumerable (r.e.) if there is a (partial or total) recursive function 
(/? : N ^ N enumerating E, that is E = {^{n) : n e N}. li E ^ 4> can be 
effectively converted into a total recursive function ip which enumerates the same 
set E. We recall a useful characterization of r.e. sets: a set C N is said to be 
semi-decidable if there is a recursive function : N ^ N whose domain is E, that 
is ^p(n) halts if and only \in & E. A set is r.e. if and only if it is semi-decidable, and 
the corresponding recursive functions can be effectively converted one another. We 
will freely use this equivalence, using in each particular situation the most adapted 
characterization. 

2.1. Algorithms and uniform algorithms. Strictly speaking, recursive func- 
tions only work on natural numbers, but this can be extended to the objects 
(thought as "finite" objects) of any countable set, once a numbering of its elements 
has been chosen. We will use the word algorithm instead of recursive function when 
the inputs or outputs are interpreted as finite objects. The operative power of al- 
gorithms on the objects of such a numbered set obviously depends on what can be 
effectively recovered from their numbers. 

More precisely, let X and Y be such numbered sets such that the numbering of 
X is injective (it is then a bijection between N and X). Then any recursive function 
(y9 : N N induces an algorithm A : X ^ Y . The particular case AT = N will be 
much used. 

For instance, the set Q of rational numbers can be injectively numbered Q = 
{goi9ij ■ • •} in an effective way: the number i of a rational a/b can be computed 
from a and 6, and vice versa. We fix such a numbering: from now and beyond qi 
will designate the rational number which has number i. 

Now, let us consider a computability notion in the real number set, here for a 
number to be computable means that there is an algorithm which can approximate 
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the number up to any error. We remark that this notion was already introduced 
by Turing in [Tur36j . 

Let a; be a real number and define Q^(a;) :— {i G N : qi < x}. 

Definition 1. We say that: 

9 X is lower semi-computable if the set Q'^{x) is r.e. 

• X is upper semi-computable if the set Q^(— x) is r.e. 

• X is computable if it is lower and upper semi-computable. 

Equivalently, a real number is computable if and only if there exists an algorith- 
mic enumeration of a sequence of rational numbers converging exponentially fast 
to X. That is: 

Proposition 1. A real number is computable if there is an algorithm ^ : N ^ Q 
such that \A{n) — x\ < 2^" for all n. 

Uniformity. 

Algorithms can be used to define computability notions on many classes of math- 
ematical objects. The definition of computability notions will be particular to a 
class of objects, but they will always follow the following scheme: 

An object O is computable if there is an 
algorithm A Y which computes O in some way. 

Each computability notion comes with a uniform version. Let (Oi)igN be a 
sequence of computable objects: 

Oi is computable uniformly in i if there is an algorithm 
: N X X — > F such that for all i, Ai = A{i, .) : X ^ Y computes Oi. 

For instance, the elements of a sequence of real numbers {xi)i^-fj are uniformly 
computable if there is a algorithm ^ : N x N ^ Q such that \A{i, n) ~ Xi\ < 2^" 
for all j, n. 

In each particular case, the computability notion may take a particular name: 
computable, constructive, eff'ective, r.e., etc... so the term "computable" used above 
shall be replaced. 

2.2. Computable metric spaces. A computable metric space is a metric space 
with an additional structure allowing to interpret input and output of algorithms as 
points of the metric space (for an introduction to other approaches to this concept 
see |HR07'). This is done in the following way: there is a dense subset (called 
ideal points) such that each point of the set is identified with a natural number. 
The choice of this set is compatible with the metric, in the sense that the distance 
between two such points is computable up to any precision by an algorithm getting 
the names of the points as input. Using this simple assumptions many constructions 
on metric spaces can be implemented by algorithms. 

Definition 2. A computable metric space (CMS) is a triple X = (X,d,S), 
where 

(i) {X, d) is a separable metric space. 

(ii) S = {sijiGN is a countable set of elements from X (called ideal points^ 
which is dense in {X,d). 

(iii) The distances between ideal points d{si,Sj) are all computable, uniformly in 
i,j (there is an algorithm that gets the names of two points and an allowed 
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error as an input and outputs the distance between two points up to the 
given approximation). 

5* is a numbered set, and the information that can be recovered from the rmmbers 
of ideal points is their mutual distances. Without loss of generality, we will suppose 
the numbering of S to be injective: it can always be made injective in an effective 
way. 

We say that in a metric space {X,d), a sequence of points {xn)neN converges 
fast to a point x if d{xn,x) < 2~" for all n. 

Definition 3. A point x G X is said to be computable if there is an algorithm 
A:N^S such that {A{n))nen converges fast to x. 

Wc define the set of ideal balls to be B := {B{si,qj) : Si G S,qj G Q>o}- We 
fix a numbering B = {Bq, Bi, . . .} which makes the number of a ball effectively 
computable from its center and radius and vice versa (this numbering may not be 
injective). B is a countable basis of the topology. 

Definition 4 (Constructive open sets). We say that an open set U is constructive 

if there is an algorithm A : N ^ B such that U = 1J„ A{n). 

Observe that an algorithm which diverges on each input n enumerates the empty 
set, which is then a constructive open set. Sequences of uniformly constructive open 
sets are naturally defined. 

Example 1. We give some example of constructive open sets: 

• The whole space X is constructive open. 

• Every finite union or intersection of ideal balls {Bn^, . . . ,Bn^} is a con- 
structive open set, uniformly in (m, . . . ,nfc). 

• // (i7i)i£N is a sequence of uniformly constructive open sets, then \J - Ui is 
a constructive open set. 

Remark 1. IfUis constructively open, belonging to U for an ideal point is semi- 

decidable: there is an algorithm A : S ^ N which halts only on ideal points belonging 
to U. Equivalently, the set of ideal points lying in U is r.e. (as a subset of S): there 
is an algorithm A •.'^ ^ S enumerating Sf\U. Hence {U, Sf\U,d) has a natural 
structure of computable metric space. 

Definition 5 (Constructive G^-set). A constructive G^-set is an intersection 

of a sequence of uniformly constructive open sets. 

Obviously, an intersection of uniformly constructive G^-sets is also a constructive 
G^-set. 

Let [X.Sx = {•sf" . s^, ...}, dx) and (Y, Sy = {sY , ■■■}:dY) be computable 
metric spaces. Let also Bf and Bf be enumerations of the ideal balls in X and 
Y. A computable function X — > F is a function whose behavior can be computed 
by an algorithm up to any precision. For this it is suflacient that the preimage of 
each ideal ball is calculated with any precision. 

Definition 6 (Computable Functions). A function T : X ^ Y is computable if 
T~^{BY) is a constructive open set, uniformly in i. That is, there is an algorithm 
A:NxN^B^ such that T-\BY) = U„ Ai, n) for all i. 

A function T : X Y is computable on D C X if there are uniformly 
constructive open sets Ui such that T~^{BY) (1 D = Uir\ D. 
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Remark 2. We remark that if T is computable then all T{sf-) are computable 
uniformly in i: there is an algorithm : N x N ^ such that 
converges fast to T{sf) for all i. 

The algorithm just semi-decides for each ideal ball in Y if sf- is contained in 
its preimage. The process will stop for each ideal hall that contains T{sf), which 
allows to extract a sequence of ideal points of Y which converges fast to T{sf). 

We remark that there are many more or less equivalent characterizations of 
computable functions between CMS, see pR07j for more details. 

The following ( |GalOO| . Def. 21 ) is a criteria to check computability of a large 
class of uniformly continuous functions. 

Proposition 2. IfT satisfies the following. ■ 

• all T{sf) are computable points, uniformly in i, 

• T is recursively uniformly continuous: there is an algorithm A : Q>o — > 
Q>o such that for all e G Q>o, d{x, x') < A{e) d{T{x), T{x')) < e, 

then T is computable. 

Proof Let E = {{i,j) G : d{T{sf),s) + qj < r}: this is a r.e. subset of N 
(uniformly in s,r) by the first condition. Then one can show that T^^{B(s,r)) — 
U(^^,)eEB{s^,A{q,)). □ 

2.3. Computable measures. When X is a computable metric space, the space of 
probability measures over X, denoted by M{X), can be endowed with a structure 
of computable metric space. Then a computable measure can be defined as a 
computable point in A4{X). 

Let X = {X, d, S) be a computable metric space. Let us consider the space 
A4{X) of measures over X endowed with weak topology, that is: 

(2.1) fj,^ —^ fj, iS fi^^f — > /i/ for all real continuous bounded / 

where fif stands for J f dfi. 

If X is separable and complete, then Ai{X) is separable and complete. Let 
D C A4{X) be the set of those probability measures that are concentrated in finitely 
many points of S and assign rational values to them. It can be shown that this is a 
dense subset (^ |Bil68] ). Let Aini,..,Tifc,mi,...mfc denote the measure concentrated over 
the finite set {sm, ■ ■ • , Sn^} with g^. the weight of Sm- 

We consider Prokhorov metric p on A4{X) defined by: 

p{fi, v) := inf{e G M+ : p{A) < J^(A') + e for every Borel set A}. 
where A'^ — {x : d{x,A) < e}. 

This metric induces the weak topology on M{X). Furthermore, it can be shown 
that the triple {Ai{X), D, p) is a computable metric space (see |Gac05j . jHR07j ). 

Definition 7. A measure p, is computable if there is an algorithmic enumeration 
of a fast sequence of ideal measures {Pn)neN C D converging to p in the Prokhorov 
metric and hence, in the weak topology. 

We need a criteria to check that a measure is computable. Let us then introduce 
(following |Gac05j ) a certain fixed, enumerated sequence of Lipschitz functions. Let 
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gs,r,e = |1- \d{x,s) - r\ + /e\ 



+ 



where s € S', r, e e Q and |a|+ — max{a, 0}. 

These are Lipschitz functions equal to 1 in the ball B{s,r), to outside B{s,r+e) 
and with intermediate values in between. It is easy to see that the real valued 
functions gsi,rj,ek '■ X —^ M. are computable, uniformly in i,j, k. 

Let be the smallest set of functions containing J^q and the constant 1, and 
closed under max, min and rational linear combinations. Clearly, this is also a 
uniform family of computable functions. We fix some enumeration vjr oi and we 
write gn for yj={n) G T . We remark that this set is dense in the set of continuous 
functions with compact support. 

The following lemma, proved in [Gac05| , shows that the approach to define com- 
putable measures we adopted, approximating measures with measures supported 
on finite ideal sets is compatible with viewing the space of measures as the dual of 
continuous functions, i.e. a measure is computable if and only if it is a computable 
function : R. 

Lemma 1. Let T ~ {51,. 92, •■•} &e the set introduced above. A probability measure 
fi is computable if and only if Jg^d/i is computable uniformly in i. 

Together with the previous lemma, the following result (see [HR07| ) will be all 
we use about computable measures: 

Lemma 2. A probability measure /i is computable if and only if the measure of 
finite union of ideal balls fJ-{Bi-^ U . . . U Bi^,) is lower semi-computable, uniformly in 
ii, . . . ,ifc. 

2.4. Computable probability spaces. To obtain computability results on dy- 
namical systems, it seems obvious that some computability conditions must be 
required on the system. But the "good" conditions, if any, are not obvious to 
specify. 

A computable function defined on the whole space is necessarily continuous. 
But a transformation or an observable need not be continuous at every point, as 
many interesting examples prove (piecewise-defined transformations, characteristic 
functions of measurable sets,...), so the requirement of being computable everywhere 
is too strong. In a measure-theoretical setting, the natural weaker condition is to 
require the function to be computable almost everywhere. It is not sufficient, and 
a computable condition on the set on which the function is computable is needed 
to assure the existence of computable points inside it. A condition which makes 
things work is the following one: 

Definition 8. A computable probability space is a pair {X, jj.) where X is a 
computable metric space and /i a computable Borel probability measure on X. 

Let Y be a computable metric space. A function (X, ^) —t Y is 5 a. e. -computable 
if it is computable on a constructive Gs-set of measure one, denoted by dom/ and 
called the domain of computability of /. 

A morphism of computable probability spaces f : {X, fi) (Y, h') is a morphism 
of probability spaces which is S a. e. -computable. 

Remark 3. A sequence of functions fn is uniformly 6 a. e.- computable if the func- 
tions are uniformly computable on their respective domains, which are uniformly 
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constructive Gs-sets. Remark that in this case intersecting all the domains provides 
a constructive Gs-set on which all fn are computable. In the following we will apply 
this principle to the iterates fn = T"' of a 5 a. e. -computable function T : X ^ X, 
which are uniformly 6 a. e. -computable. 

3. Constructive Borel-Cantelli sets 

Given a space X endowed with a probability measure fi, the well known Borel 
Cantelli lemma states that if a sequence of sets Ak is such that J2 M(^fc) < c» then 
the set of points which belong to finitely many A/.' a has full measure. In this section 
we show that if the A^ are given in some "constructive" way (and n is computable) 
then this full measure set contains some computable points. Hence this set contains 
points which can effectively be constructed. 

Definition 9. A sequence of positive numbers ai is effectively summable if the 
sequence of partial sums converges effectively: there is an algorithm A : Q — > N 
such that if A{e) = n then J2i>n ^« — ^• 

Remeirk 4. A computable sequence of positive real numbers is effectively summable 
if and only if its sum is a computable real number. 

For sake of simplicity, we will focus on the complements Un of the A„. 

Definition 10. A constructive Borel-Cantelli sequence is a sequence ([/n)rteN of 
uniformly constructive open sets such that the sequence (/u(f/^)) is effectively sum- 
mable. 

The corresponding constructive Borel-Cantelli set is r\n>k 

The Borel-Cantclli lemma says that every Borel-Cantelli set has full-measure: we 
are going to see that every constructive Borel-Cantelli set contains a dense subset 
made of computable points. 

Lemma 3 (Normal form lemma). Every constructive Borel-Cantelli sequence can 
be effectively transformed into a constructive Borel-Cantelli sequence (?7„)neN giving 
the same Borel-Cantelli set, with tiiU^) < 2~". 

Proof, consider a constructive Borel-Cantelli sequence (V^). As /u(X \ V^) is effec- 
tively summable, an increasing sequence (ni)j>o of integers can be computed such 
that for all i, Y.n>n, \ K) < 2-\ 

We now gather the Vn by blocks, setting: 

Ui := fl Vn 
Ui is constructively open uniformly in i, and: 

M(c/f ) < 2-' and u n ^» = u n = u n 

k n>k i n>ni i j>i 

□ 

In the sequel we will always suppose that a constructive Borel-Cantelli sequence 

is put in this normal form. 

Proposition 3. Every finite intersection of constructive Borel-Cantelli sets is a 
constructive Borel-Cantelli set. 
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Proof, let ([/„) and {Vn) be two constructive Borel-Cantelli sequences in normal 
form. It is easy to see that: 

u n ^ u n ^» = u n ^ k 

k n>k k n>k k n>k 

and fJ.{{Un n Vn)'^) < 2^"+^ which is effectively summable. □ 

As every effectivity notion, the notion of constructive Borel-Cantelli set naturally 
comes with its uniform version. 

Proposition 4. The intersection of any uniform family of constructive Borel- 
Cantelli sets contains a constructive Borel-Cantelli set. 

Proof, suppose that Ri — IJj, nn>fc normal form. Consider a simple bijec- 

tion : {(n, i) : < « < n} ^ N (for instance, (p{n,i) = n{n+l)/2 + i) computable 
in the two ways and define the sequences (V^m)meN and {a,in)rnefi by Vm = and 
a„i ~ 2~" where i) = m. Obviously > a.m. 

A simple calculation shows that "^m — 4. So (V^m) is a constructive Borel- 
Cantelli sequence. 

Fix some i. If a point is outside J/^ for infinitely many n, it is outside Vm for 
infinitely many m. That is to say: 

u n K.c|jn^^-^^ 

k m>k k n>k 

As it is true for every i, the constructive Borel-Cantelli set induced by {Vm)m is 
included in every Ri. □ 

3.1. Computable points in constructive Borel-Cantelli sets. The Borel- 
Cantelli lemma can be strengthened for constructive Borel-Cantelli sequences ob- 
taining that they contain computable points. 

Theorem 1. Let X be a complete CMS and fi a computable Borel probability 
measure on X . 

For every constructive Borel-Cantelli set R, the set of computable points lying 
in R is dense in the support of fi. 

In order to the prove this theorem, we need the following lemma to construct a 
computable point from what could be called a shrinking sequence of constructive 
open sets. 

Lemma 4 (Shrinking sequence). Let X be a complete CMS. Let Vi be a sequence 
of non-empty uniformly constructive open sets such that Vi+i C Vi and diam(Vi) 
converges effectively to 0. Then Hi ^ singleton containing a computable point. 

Proof. As Vi is non-empty there is a computable sequence of ideal points Si € Vi. 
This is a Cauchy sequence, which converges by completeness. Let x be its limit: it 
is a computable point as diam(V^;) converges to in an effective way. Fix some i: 
for all j > i, Sj &Vj(=ViSOx — hmj^oo sj e Vi. Hence x e f^^Vi = f]^Vi. □ 

Proof of theorem]^ Let {Un)n be a constructive Borel-Cantelli sequence, in normal 
form {n{Un) > 1 — 2~", see lemma [3]). Let B be an ideal ball of radius r < 1 and 
positive measure. In B we construct a computable point which lies in [j^ r\k>n ^k- 
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To do this, let Vq = B and no be such that fJ.{B) > 2^""+^ (such an no can be 
effectively found from B): from this we construct a sequence {Vi)i of uniformly con- 
structive open sets and a computable increasing sequence {ni)i of natural numbers 
satisfying: 

(1) m) + f^ink>nUk) > 1, 

(2) V, C f]n„<k<n. Uk, 

(3) diam(y,) < 2-*+i, 

(4) F,+i c V,. 

The last two conditions assure that Hi is a computable point (lemma |4]) , the 
second condition assures that this point lies in nfe>ri Uk- 
Suppose Vi and rij have been constructed. 

Claim 1. There exist m > ni and an ideal ball B' of radius 2~*~^ such that 

(3.1) fi{v,n fl [/fens') > 2-™+!. 

We now prove of the claim: By the first condition, ii{Vi np|j,>^. Uk) > so there 
exists an ideal baU B' of radius 2"'"i such that /i(Vj n C\k>ni Uk n B') > 0. There 
is m > ni such that /i(Vi n nfe>n Uk n B') > 2~™+^, which impHes the assertion, 
and the claim is proved. 

As inequality p.ip can be semi-decided, such an m and a i?' can be effectively 
found. For T^i+i, take any finite union of balls whose closure is contained in Vi D 
r\n <k<mUk H B' and whose measure is greater than 2^™+^. Put rti+i = to. 
Conditions 2., 3. and 4. directly follow from the construction, condition 1. follows 
from /i(Vi+i) > 2"™+-'^ > 1 — A'(nfe>m Uk) (the sequence is in normal form). □ 

The following corollary allows to apply the above criteria to a uniform infinite 
sequence of constructive Borel-Cantelli sets. 

Corollary 1. Let X be a complete CMS and fi a computable Borel probability 
measure on X . 

Given a uniform family {Ri)i of constructive Borel-Cantelli sets, the set of com- 
putable points lying in p|,- Ri is dense in the support of fi. 

Proof, this a direct consequence of proposition 3] and theorem [T] □ 

We remark that, in the particular case of Cantor spaces with an uniform mea- 
sure a result of this kind can also be obtained from |Sch71| since it is possible 
to relate Borel Cantelli sequences to Schnorr tests. This relation is developped in 
[GHROSj giving some new connections between Schnorr randomness and dynamical 
typicality. 

3.1.1. Application to convergence of random variables. Here, {X, fi) is a computable 
probability space. 

Definition 11. A random variable on {X^n) is a measurable function / : X — > 
R. 

Definition 12. Random variables fn effectively converge in probability to f 

if for each e > 0, ^{x : \fnix) — f{x)\ < e} converges effectively to 1, uniformly 
in e. That is, there is a computable function n{t,5) such that for all n > n{e,6), 
-f\>e]<S. 
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Definition 13. Random variables /„ effectively converge almost surely to f 

i!/supj.>„ l/ri — /I effectively converge in probability to 0. 

Theorem 2. Let /«,/ be uniformly S a. e.- computable random variables. If fn 
effectively converges almost surely to f then the set {x : fn{x) —* f{x)} contains a 
constructive Borel-Cantelli set. 

In particular, the set of computable points for which the convergence holds is 
dense in Supp(/x). 

Proof. Let D — p|„ Z?„ be a constructive G^-set of full measure on which all /„, / 
are computable. _D„ are uniformly constructive open sets, and we can suppose 
Dn+i ^ Dn (otherwise, replace -D„ by Dq n . . . n I?„). 

There are uniformly constructive open sets t/„(e) such that [/„(e) D — [|/n — 
/I < e] n D. y^(nn>fc Un{e)) converges effectively to 1, uniformly in e so it is 
possible to compute a sequence (fc^); such that Ai(n„>fc C^n(2~*)) > 1 — 2~* for all 
i. Put Vi = H/c <n<fe +1 tAi(2^*) n Di'. Vi is constructively open uniformly in i and 
niYi) > 1 — 2~*. The sets Vi form a constructive Borel-Cantelli sequence, and if a 
point X is in the corresponding Borel-Cantelli set then x Cz D and there is zq such 
that X £ Vi for all i > io, so \fn{x) ~ f{x)\ < 2^* for all n > ki,i > io. Hence 



Let X be a metric space, let T : X i-^ X be a Borel map. Let /i be an invariant 
Borel measure on X, that is: IJ,{A) = iJ,{T~^{A)) holds for each measurable set A. 
A set a A is called T-invariant if T~^{A) = A(mod 0). The system {T,iJ,) is said 
to be ergodic if each T-invariant set has total or null measure. In such systems 
the famous Birkhoff ergodic theorem says that time averages computed along fj, 
typical orbits coincides with space average with respect to /n. More precisely, for 
any / G L^{X) and it holds 



for /i almost each x, where Sl = f + foT+... + fo T"^^. 

If a point X satisfies equation 14. II for a certain /, then we say that x is typical 
with respect to the observable /. 

Definition 14. // x is typical w.r.t any continuous function / : X — > R with 
compact support, then we call it a /x- typical point. 

In this section we will see how the constructive Borel-Cantelli lemma can be used 
to prove that in a large class of interesting systems there exists computable typical 
points. 

Let us call {X, /i, T) a computable ergodic system if {X, ^) is a computable 
probability space, T is an endomorphism (i.e. an 6 a. e. -computable measure- 
preserving transformation) and (X, fi, T) is ergodic. 

Before to enter in the main theme of typical statistical behaviors let us see 
an easier topological result in this line. One of the features of undecomposable 
(topologically transitive) chaotic systems is that there are many dense orbits, the 
following shows that if the system is computable then there are computable dense 
orbits. 



fn{x) f{x). 



□ 



4. Pseudorandom points and dynamical systems 



(4.1) 
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We remark that this result can also be obtained as a corollary of the constructive 
Baire theorem |YMT99) . 

Theorem 3. Let X be a computable complete metric space and T : X ^ X a 
transformation which is computable on a dense constructive open set. If T has a 
dense orbit, then it has a computable one which is dense. 

In other words, there is a computable point x € X whose orbit is dense in X. 
Actually, the proof is an algorithm which takes an ideal ball as input and computes 
a transitive point lying in this ball. 

Proof. {Bi)i^tii being an enumeration of all ideal balls, define the open sets Ui = 
dom(f) n |J„ T~^Bi which are constructive uniformly in i. By hypothesis, Ui is 
also dense. Ui is the set of transitive points. From any ideal ball B{so,ro) we 
effectively construct a computable point in i3(so, tq) n - Ui. 

If B{si,ri) has been constructed, as Ui is dense B{si,ri) O Ui is a non-empty 
constructive open set, so an ideal ball B{s,r) C B{si,ri) n Ui can be effectively 
found (any of them can be chosen, for instance the first coming in the enumeration). 
We then set i?(si+i, ri+i) := B{s,r/2). 

The sequence of balls computed satisfies: 

B{si+i,r,+i) C B{s,,ri) n Uq H . . . H U, 

As (r'i)igp} is a decreasing computable sequence converging to and the space 
is complete, (si)igN converges effectively to a computable point x. Then {x} = 
aBis,,r,)Cf]^U,. □ 

4.1. Computable typical points. We will use the results from the previous sec- 
tion to prove that computable typical points exist for a class of dynamical systems. 
Each time the set of typical points is a constructive Borel-Cantelli set, theorem [5] 
applies. 

For instance, in the case of the shift on the Cantor space with a Bernoulli mea- 
sure, the Birkhoff ergodic theorem reduces to the strong law of large numbers, 
which proof is simpler and makes explicit use of the Borel-Canteli lemma. This is 
possible thanks to the independence between the random variables involved, but 
strict independence is actually unnecessary: the proof can be adapted whenever 
the correlations between the random variables decrease sufficiently fast. 

Definition 15. We say that a .system {X,T,fj,) is In^-ergodic for observables in 



some 



set of functions B if for each [ipTip) S there is c^,^ > such that 



^ -<n J J 



(In(nj)^ 



Now we can state: 

Theorem 4. Let {X,T,fj,) be a dynamical system which is hi^ -ergodic for observ- 
ables in some set B of bounded observables. For each <j) £ B, the almost-sure 
convergence: 

is effective. 



CONSTRUCTING PSEUDO-RANDOM POINTS 



15 



Note that for the moment, no computabihty assumption is needed on the system. 
As announced, the proof is an adaptation of the proof of the strong law of large 
numbers. We first prove two lemmas. 

Lemma 5. There exists a computable sequence rii such that: 

• /9j := converge effectively to 1, 

• Y^j^^^yi is effectively summahle. 

Proof. For instance, take Ui = [(1 + with < a < 1/2. 

From now on, we denote — by /„. 



□ 



Lemma 6. The almost-sure convergence of the subsequence fm to J(f>dfi is effec- 
tive. 



Proof. For 5 > 0, define the deviation sets: 



fn{x) - (("in 



> S 



xeX : 
By Tchebytchev inequality, 

S'l^iAtiS)) < 



Since adding a constant to cj) does not change this quantity, without loss of 
generality, let us suppose that /0d/i = 0. Then 



/„ - (j) dfi 



L2 



^0 



d/i 



■ oT+... + 0oT 



ji-i 



d/^ 



by invariance of /x this is equal to 

^2 E <t>-T'-'4>)d^^ 



i<j<n 



hence. 



As ^ 



< 



n 
n 



n ^ 

ln(n)2 



>oT''(t>dn 



In(ni)' 



is effectively summable (by choice of n^, see lemma [5]) uniformly 



in (5, it follows that /„; converge effectively almost-surely to Jcjidfi. 



□ 



As n.i is not dispersed too much, the almost-sure convergence of the subsequence 
/„. implies that of the whole sequence /„. Actually the effectivity is also preserved. 
We now make this precise. 

Lemma 7. For Ui < n < Ui+i and Pi :— one has: 

(4.2) ll/n.-/n|loo<2(l-/?,)|l0|U. 
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Proof. Let M = To see this, for any k,l,P with (3 < k/l < 1: 

St_Sl_ ^ f k\St St,oT^~' 
k I \ I J k I 

< (1 - /3)M + ^ = 2(1 - /3)M, 

Taking (3 = (3^ and k = rii, I = n first and then k = n, I = n.i^i gives the result. □ 

Proof of theorem^ Let 6,e > 0. To prove that /« converge effectively almost- 
surely, one has to compute some p (from 6 and e) such that Ai(Un>p < ^• 

As /3j converge effectively to f, one can compute such that if i > «o then 
/3j > f - 5/(4M). Inequality 221 then implies 

U ^„(5) C A..(<5/2). 
Indeed if < n < n^+i and |/ni(x) — J(j)dfi\ < 6/2 then |/„(a;) — J(j)dn\ < 

\fn{x) - fnM\ + IfnM " UM < 6. 

As /„. converge effectively almost-surely, one can compute some jo such that 
m(Uj>jo < Let p = Uk where k = max(io,jo): Un>pAi('5) ^ 

Uj>jo whose measure is less than e. □ 

Corollary 2. Lei {X,T,fi) be a computable dynamical system which is It? -ergodic 
for observables in some set B of bounded functions and let (p be a S a. e.- computable 
observable in B. 

The set of points which are typical w.r.t <j) contains a constructive Borel-Cantelli 
set. In particular, it contains computable points. 

Proof. Apply theorem [5] to the sequence of uniformly 5 a. e. -computable functions 
/„ — ^ which converge effectively almost-surely by theorem 21 □ 

Remark 5. In the proof of thm^ we see that the constructive Borel-Cantelli set 
depends in an effective way on {{(t^W^ and c^^^. This gives the possibility to operate 
in a way to apply Prop. and Cor. [7] to find a constructive Borel Cantelli set and 
computable points contained in the set of points typical with respect to a uniform 
family (j>„Ti. 

By the above remark, to construct /^-typical points (see definition I14p using the 
above mentioned results, the following conditions are sufficient: 

Theorem 5. If a computable system is hi' -ergodic for observables inT = {gi, 52, ■ ■ •} 
(this set was defined in section \2. S\] and the associated constants Cg. (see definition 
\16\} can be estimated uniformly in i ( there is an algorithm A : N — > Q such that 
A{i) > Cg.) then it has a set of computable ^-typical points which is dense in the 
support of pL. 

Proof. We remark that T is dense in the set of continuous functions on X with 
compact support (with the sup norm) hence a computable point which is typical 
for each gi is /i- typical. Such points can be found by applying theorem [2] for each 
gi and using proposition 2] as explained in remark [5] □ 
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4.1.1. In -mixing. We will apply this to systems having a stronger property: they 
are mixing, with logarithmical speed. More precisely, this can be quantified using 
the correlation functions: 



0oT"-0d/i— / (/)d/i /V'd/i 



which measures the dependence between observation through and -0 at times 
n ^ 1 and respectively (possibly with ip = 0). Note that Cn{(f>,ip) — corre- 
sponds, in probabilistic terms, to ^oT" and being independent random variables. 

Definition 16. We say that a system {X,T,ii) has In^-decay of correlations 

for observables in some set of functions B if for each (0, "0) G there is c^^^ > 
such that 



(ln(n))2 



for all n > 2. 



Lemma 8. If a system has In -decay of correlation for observables in B then it 
is It? -ergodic for observables in B. The ergodicity constants depend in an effective 
way on the mixing constants. 

Proof. We first prove that for all rt > 2, 



(4.3) 

For n > 56, 



E 

fe=2 



1 



< 



2n 



ln(n) 



1 



E in(fc)2 



< 



< 



x—55 



dx 
ln(x)2 



x— 55 

2n 



1 



ln(a;)2 
110 



ln(a;)-^ 



da; (as 55 > ln(4)) 



which, combined with < 10 and 



Finally, for n > 2 

n ^ — ' 



ln(n)2 ln(55)2 

5 1 
^k=2 ln(fe)2 



i^]l"yi > 6, gives inequality 



o T*0 d/i — / d/i ipdfi 



< 



< 



< 



i^a(0,0) 



z<n 

ln(n)2 



4c, 



□ 



4.2. Application: computable absolutely normal numbers. An absolutely 
normal (or just normal) number is, roughly speaking, a real number whose digits 
(in every base) show a uniform distribution, with all digits being equally likely, all 
pairs of digits equally likely, all triplets of digits equally likely, etc. 

While a general, probabilistic proof can be given that almost all numbers are 
normal, this proof is not constructive and only very few concrete numbers have 
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been shown to be normal. It is for instance widely believed that the numbers 
•\/2, TT and e are normal, but a proof remains elusive. The first example of an 
absolutely normal number was given by Sierpinski in 1916, twenty years before the 
concept of computability was formalized. Its construction is quite complicate and 
is a priori unclear whether his number is computable or not. In |BF02j a recursive 
reformulation of Sierpinski's construction (equally complicate) was given, furnishing 
a computable absolutely normal number. 

As an application of theorem[2]we give a simple proof that computable absolutely 
normal numbers are dense in [0, 1]. 

Let b be an integer > 2, and Xt the space of infinite sequences on the alphabet 
Sb = {0, . . . , 6 — 1}. Let T = cr be the shift transformation on Xb, and A be the 
uniform measure. A real number r G [0, 1] is said to be absolutely normal if for all 
& > 2, its &-ary expansion G Xf, satisfies: 

1 ""^ 1 
lim - o cr'irb) = — r for all w S S^. 

4 = 

Theorem 6. The set of computable reals which are absolutely normal is dense in 
[0,1]. 

Proof, for each base b > 2, consider the transformation Tf, : [0, 1] — > [0, 1] defined 
by Tb{x) = 6x(mod 1). The Lebesgue measure A is T;,-invariant and ergodic. The 
partition in intervals [k/b, [k + l)/b[ induces the symbolic model (S^, cr. A) which 
is measure-theoretically isomorphic to ([0, 1], T;,, A): the interval [k/b, (k + l)/b[ is 
represented by G T,b. For any word ui e E;, define I{w) to be the corresponding 
interval [O.w, O.uJ + 2" !"'!]. 

Defining domTb := [0, 1] \ : < k < b} (the interior of the partition) makes 
Tb SL 6 a. e. -computable transformation. The observable f^ :— l/(^) is also S a.e.- 
computable, with dom/tu = [0, 1] \ dl{w). 

Actually, since fw o c" and f^ are independent for n > \w\, theorem [2] applies 
to ([0, l],Th,A) and fw Therefore, the set of points (for the system (T^, A)) which 
are typical w.r.t the observable fw contains a constructive Borel-Cantelli set Rb,w 
Furthermore, Rb.w is constructive uniformly in 6, w G Sfc. Hence, by corollary [U 
their intersection, which is made of absolutely normal numbers, contains a dense 
set of computable points. □ 

5. Dynamical systems having computable typical points 

We will see that in a large class of dynamical systems which have a single physi- 
cally relevant invariant measure, the computability of this measure and related Cg- , 
for observables in can be proved, hence we can apply Thm. [5] to find pseudoran- 
dom points in such systems. 

5.1. Physical measures. In general, given {X,T) there could be infinitely many 
invariant measures (this is true even if we restrict to probability measures). Among 
this class of measures, some of them are particularly important. Suppose that we 
observe the behavior of the system {X, T) trough a class of continuous functions 
fi : X — > M. We are interested in the statistical behavior of fi along typical orbits 
of the system. Let us suppose that the time average along the orbit of x exists 

A,(/,)= hm i^/,(r"(x)) 
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this is a real number for each ft. Moreover Ax{fi) is linear and continuous with 
respect to small changes of fi in the sup norm. Then the orbit of x acts as a 
measure ^.^ and Ax{fi) = Jfid^^ (moreover this measure is also invariant for 
T). This measure is physically interesting if it is given by a "large" set of initial 
conditions. This set will be called the basin of the measure. If X is a manifold, it is 
said that an invariant measure is physical (or SRB from the names of Sinai, Ruelle 
and Bowen) if its basin has positive Lebesgue measure (see |You02j for a survey 
and more precise definitions). 

In what follows we will consider SRB measures in the classes of systems listed 
below, 

(1) The class of uniformly hyperbolic system on submanifolds of R". 

(2) The class of piecewise expanding maps on the interval. 

(3) The class of Manneville-Pomeau type maps (non uniformly expanding with 
an indifferent fixed point). 

All these systems, which are rather well understood, have a unique physical 
measure with respect to which correlations decay is at least polynomial. Further- 
more, in each case, the corresponding constants can be estimated for functions in 
J- . The computability of the physical measures is proved case by case, but it is 
always a consequence of the fact that, in one way or another, the physical measure 
is "approached" by iterates of the Lebesgue measure at a known speed. 

5.2. Uniformly hyperbolic systems. To talk about SRB measures on a system 
whose phase space is a manifold, we have to introduce the Lebesgue measure on a 
manifold and check that it is computable. 

5.2.1. Computable manifolds and the Lebesgue measure. For simplicity we will not 
consider general manifolds but submanifolds of R". 

Definition 17. Let M be a computable metric subspace o/R". We say that M 
is a m dimensional computable submanifold of R" if there exists a computable 
function f : M x B{0, 1) ^ M (where B(0, 1) is the unit ball of M."" and M x B{0, 1) 
with the euclidean distance is a CMS in a natural way) such that for each x e AI , 
fx — f{x, ■) is a diffeomorphism with all k derivatives being computable. 

For each x, the above fx is a map whose differential at any z £ B{0, 1) is a 
linear, rank m function Df^.z '■ R™ R"- This can be seen as a composition of 
two functions Dfx,z = Dfl^o Dfl ^ such that Dfl ., : R™ R™ is invertible and 
Dfl ., : R" ^ R" is an isometry. 

Let us denote the image of 5(0,1) by f^. Then the Lebesgue measure of 
D C Bx is defined as 



This does not depend on the choice of Bx and fx, and it give rise to a finite 
measure (Lebesgue measure) on M (see jGMS98j page 74). This measure is indeed 
the m dimensional Hausdorff measure on M. Moreover, the Lebesgue measure is a 
computable measure. 

Lemma 9. The Lebesgue measure on a computable submanifold ofW'' is com- 
putable. 
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Proof (sketch). Suppose that A is a constructive open subset of some Bg, where s is 
an ideal point of M. Since the function det{D f^ ^) is computable and the function 
1 f-i-(A) i^) is lower semi-computable, we can lower semi-compute the value m{A). In 
particular, there is a base of ideal balls whose measures are lower semi-computable. 
Let B and B' be such balls. Since these balls have zero measure boundaries, 
we can compute the measure of their intersection (which is a constructive open 
included in B). Hence any constructive open set can be decomposed into a (same 
measure) disjoint union of constructive open sets whose measures can be lower 
semi-computed. By lemma [21 m is computable. □ 

5.2.2. The SRB measure of uniformly hyperbolic systems. Let us consider a con- 
nected computable manifold A/. Let us consider a dynamical system (Af, T) 
where T is a computable diffeomorphism on M. 

Let us consider a constructive open forward invariant set Q G M (i.e. T{Q) C 
Q). Let us consider the (attracting) set 

A= n T"(Q). 

n>0 

Suppose that A contains a dense orbit and that it is an hyperbolic set for T, which 
means that the following conditions are satisfied. 

There is a splitting of the tangent bundle of M on A: Tt^M = Ej^® E^ (at each 
point X of A the tangent space at x can be splitted in a direct sum of two spaces, 
the stable directions and the unstable ones) and a Ao < 1 such that 

• the sphtting is compatible with T, that is: DT^iE^) = and DT-^{E^) 

• The dynamics expand exponentially fast in the unstable directions and 
contracts exponentially fast in the stable directions in an uniform way, that 
is: for each x £ K and for each v E E^ and w S -B", \DTx{v)\ < \Xov\ and 
\DT-\w)\ < \Xow\. 

Under these assumptions it is known that 

Theorem 7. (see |Via97j e.g.)There is a unique invariant SRB measure /i sup- 
ported on A. Moreover the measure is ergodic and its basin has full Lebesgue mea- 
sure on Q. 

This measure has many good properties: it has exponential decay of correlations 
and it is stable under perturbations of T (see |Via97| e.g.). Another good property 
of this measure is that it is computable. 

Theorem 8. If M and T are , computable and uniformly hyperbolic as above, 
then the SRB measure fi is computable. 

Proof. Let m be the Lebesgue measure on Q normalized by m{Q) = 1, clearly it is 
a computable measure. From |Via97j (P rop. 4.9, Remark 4.2) it holds that there 
are A < 1 such that for each i^- Holder {v e (0, 1]) continuous observable it holds 



-0 o T" dm — ip dfi 



< A' 



where = C /IV'I dm -I- HV'lljyj where C is independent from tjj and then can be 
estimated for each uniform sequence tp^ E T uniformly in i. This means that 
for each ip^^ E J- its integral with respect to /i can be calculated up to any given 
accuracy, uniformly in i. Indeed if we want to calculate Jip^ d/i up to an error of 
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e we calculate . up to an error of e (this error is not really important as we will 
see immediately) and choose an n such that c^.A" < |. 

By this we know that ° 2^" dm — /-(/'i d/i| < |. Now we have to calculate 
JV'i ° dm up to an error of | and this will be the output. By lemma [T] then /j, 
is computable. □ 

Corollary 3. In an unif. hyp. computable system equipped with its SRB measure 
as above, the set of computable ^—typical points is dense in the support of ^. 

Proof, fi is computable by the previous theorem, and the correlations decay is given 
by proposition 4.9 in |Via97j from which follows that there is A < 1 such that for 
each {gi,gj) G J-^ it holds, 



Jg,o T^-gj dm- j gi dfi J g^ d^ 



where Cg._gj = C{J\gi\ dm + |i.gi|l | dm + ||.gj||]^) (C is a constant independent 
of gi € 11*11]^ is the Lipschitz norm, since functions in J- are Lipschitz) are 
computable uniformly in Then the result follows from theorem [S] □ 

5.3. Piecewise expanding maps. We introduce a class of discontinuous maps on 
the interval having an absolutely continuous SRB invariant measure. The density of 
this measure has also bounded variation. We will show that this invariant measure 
is computable. 

Let / be the unit interval. Let T : / / we say that T is piecewise expanding 
if there is a finite partition P — {/i, Ik} of /, such that li are disjoint intervals 
and: 

(1) the restriction of T to each interval can be extended to a monotonic 
map defined on li and the function h : I ^ R defined by h{x) — |DT(a;)|~"'^ 
has bounded variation. 

(2) there are constants C > and a > 1 such that \DT^{x)\ > Ccr" for every 
n > 1 and every x I for which the derivative is defined. 

(3) For each interval J C I there is n > 1 such that f'^{J) — I. 

We remark that by point 1), in each interval li the map is Lipschitz. We remark 
that this restriction is not strictly necessary for what follows (see }GHR| ) . we sup- 
pose it for the seek of simplicity. As said before, by classical results this kind of 
map has an absolutely continuous invariant measure (see [Via97j . chapter 3 e.g.). 

Theorem 9. IfTa piecewise expanding map as above, then it has a unique ergodic 
absolutely continuous invariant measure (x. The basin of this SRB measure has full 
Lebesgue measure. Moreover /i can be written as fi — (pm where 4> has bounded 
variation and m is the Lebesgue measure. 

Moreover as before, the SRB measure is also computable 

Proposition 5. IfT is a m— computable piecewise expanding map satisfying points 
1 ),..., 3) above then its SRB measure is computable. 

Proof. Let us consider ip Q J- and Jip o T" dm. Since T is ^a.e computable for 
the Lesbegue measure, the ends of the intervals Ii,...,Ik are computable, and so 
is (uniformly on n) for T". If T is Z-Lipschitz in each interval li then T" is Z"- 
Lipschitz in each of its continuity intervals and the Lipschitz constant of V' ° T" 
on each continuity interval can be estimated. By this is straightforward to show 
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that also Jip o T" dm can be calculated up to any accuracy (approximate Tp o T" 
by piecewise constant functions and estimate error by Lipschitz constant). 

Now, from [Via97j proposition 3.8, remark 3.2 it holds that there are A < 1, C > 
such that for each ip G 

y V o r" dm - y 

This means that since J^p o T" dm can be calculated up to any given error then 
the integral Jip d/i with respect to /i can be calculated up to any given accuracy 
and by lemma [T] then fi is computable. □ 

As Unif. Hyperbolic systems, also Piecewise Expanding maps can be shown 
to have exponential decay of correlations on bounded variation observables (see 
|Via97| Remark 3.2) and BV norm of functions in can be estimated. Hence as 
in the previous section we obtain 

Corollary 4. In an ni- computable piecewise expanding system equipped with its 
SRB measure, the set of computable typical points is dense in [0, 1]. 

5.4. Manneville-Pomeau type maps. We say that a map T : [0, 1] [0, 1] is a 

Manneville-Pomeau type map ( MP map ) with exponent s if it satisfies the following 
conditions: 

(1) there is c G (0, 1) such that, if /q = [0, c] and Ii = (c, 1], then T\f^^ ^-j and 
T|^^ extend to diffeomorphisms, which is for a; > 0, T(/o) = [0, 1], 
T(/i) = (0, 1] and T(0) = 0; 

(2) there is A > 1 such that T' > A on /i, whereas T' > 1 on (0, c] and 
r(0) = 1; 

(3) the map T has the following behaviour when x 0^ 

T{x) = X + rx^+'' {1 + u{x)) 

for some constant r > and s > and u satisfies u(0) = and u'{x) = 
0{x*~^) as a; ^ 0+ for some t > 0. 

In [Iso03j (see also jGou04| ) it is proved that for < s < 1 these systems have a 
unique absolutely continuous invariant measure, whose density / is locally Lipschitz 
in a neighborhood of each x > (the density diverges at x = 0) the system has 
polynomial decay of correlations for (1 — s)-H61der observables. Moreover we have 
that: 

Theorem 10. If T is a computable MP map then its absolutely continuous invari- 
ant measure is computable. 

Proof. Let / be the density of /i. T is topologically conjugated to the doubling 
map X ^ 2x (mod 1) hence for each small interval / there is fc > such that 
T^{I) = [0,1]. Since / is locally Lipschitz, there is a small interval J on which 
/ > i5i > 0. Let n be such that T"'{J) = [0,1]. Let / be some small interval, 
then there exist J' C J such that T'^{J') — I. Since T is A-Lipschitz, we have 
m(J') > By this, n{J') > and by the invariance of ^, ^ > ^ and 

then, as / is arbitrary, for each x € [0, 1] we have f{x) > ji- > 0. In particular, 
J is (1 — s)-H61der. Now we use the fact that the system has polynomial decay of 
correlations for (1 — s)-H61der observables. Let us consider e then we have 



<CA" ll^ll^.. 
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that jdfi = dm and J j d/i = 1, hence, by the decay of correlation of this kind of 



maps 



y"(/) o T" dm - y^d^ = J(j)o T'j d^i- Jcl^dn 



<C\\4>\W_ 



The norm 



can be estimated for functions in T , and then, as in the previous 



examples we have a way to calculate J (p dfi for each e J-" and again by lemma [U 
is computable. □ 

Corollary 5. In a computable Manneville-Pomeau type system, the set of com- 
putable typical points is dense in [0, 1]. 
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